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I. INTRODUCTION 

A fundamental challenge of high energy particle 
physics is to understand the spin structure of protons, 
neutrons, and nuclei in terms of their parton constituents. 
The increasing precision of experimental data on inclu- 
sive polarized deeply inelastic scattering (DIS) of leptons 
from nucleons allows us to perform incisive QCD analy- 
ses of polarized structure functions to reveal the spin de- 
pendent partonic structure function of the nucleon. Po- 
larized DIS lepton-nucleon scattering experiments have 
been performed at CERN, SLAC, DESY and JLAB [H- 
[l3| , and these processes have played a key role in our 
understanding of QCD and the spin structure of the nu- 
cleon There are several comprehensive analyses 
of the polarized DIS data in the literature P^-HH; this 
work provides a detailed picture of the spin structure of 
the nucleons. 

The new precision experimental data from the HER- 
MES and COMPASS collaborations [H El of the spin 
structure function g\ provides additional information 
that we shall use to study the spin structure and quark 
helicity distributions. We shall choose an approach based 
on the expansion of orthogonal polynomials; specifically, 
we will implement Jacobi polynomials as we use exper- 
imental data for each bin of Q 2 separately [43[. Previ- 
ously we applied the Jacobi polynomials to deter- 
mine the polarized valon distributions using only the pro- 
ton experimental data. In this analysis, both the unpo- 
larized and polarized valon distributions were extracted, 
so more unknown parameters were required as compared 
to the present analysis. The Jacobi polynomial expan- 
sion has also been applied to a variety of QCD analyses 
HB-liil, including the case of polarized PDFs [II EH73| . 

In the present study, we perform a NLO QCD analysis 
of the polarized deep-inelastic data in the MS- 

scheme and extract parameterizations of the polarized 
PDFs and structure functions. In Section II, we provide 
an overview of the Jacobi polynomials approach. In Sec- 
tion III we review the parametrization and evolution of 
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the PDFs. In Section IV we present the results of our fit 
to the data, and in Section V we compute the associated 
structure functions and sum rules. Section VI contains 
the conclusions. We also provide an Appendix which de- 
scribes the FORTRAN-code which is available. 



II. THE JACOBI POLYNOMIAL METHOD 

We perform a NLO fit of the polarized parton distribu- 
tions (PPFDs) using Jacobi polynomials to reconstruct 
the x dependent quantities from their Mellin moments. 
The use of Jacobi polynomials has a number of advan- 
tages; specifically, it will allow us to factorize the x and 
Q 2 dependence in a manner that allows an efficient pa- 
rameterization and evolution of the structure functions. 

For example, if we consider the spin structure function 
xgi(x, Q 2 ), we can expand this as: 

N max 

x 9l (x,Q 2 ) = x^l~xr a n (Q 2 )Q^(x) . (1) 

n=0 

Here, <d"'^(x) are Jacobi polynomials of order n, and 
N max is the maximum order of our expansion. In this in- 
stance, the Jacobi polynomials allow us to factor out the 
essential part of the ^-dependence of the structure func- 
tion into a weight function [4^|, and the Q 2 -dependence 
is contained in the Jacobi moments a n {Q 2 ). 

To be more specific, the x-dependence of the Jacobi 
polynomials can be written as 

n 

e^(*) = 5>fW)^, (2) 

where the (a, /3) coefficients are combinations of T- 
functions involving {n, a, /?}. The Jacobi polynomials 
satisfy an orthogonality relation with weight function 
x@(l — x) a as follows: 

[ dx x?(l - x) a Q^(x)Q^(x) = S k , . (3) 
Jo 

Thus, given the Jacobi moments a n {Q 2 ), the polarized 
structure function xg\{x, Q 2 ) may be reconstructed from 
Eq. (TTJ) g|. 

We can compute the Jacobi moments a n (Q 2 ) using the 
orthogonality relation to invert Eq. JTJ) to obtain: 

a n (Q 2 ) = / dxx gi (x,Q 2 )e^(x) 
Jo 

n 

= ^4 n) (a,/3) M[ X9l ,j + 2] . (4) 

3=0 

In Eq. Q, we have substituted Eq. (PJ) for xgi(x, Q 2 ) and 
introduced the Mellin transform: 



M[x5i,iV] ee f dx x N ~ 2 x gi (x,Q 2 ) . (5) 
Jo 

We can now relate the polarized structure function 
xgi(x, Q 2 ) with its moments (44| 

x 9l (x,Q 2 ) = x^(l-x) a £ 

n 

x Y, c f\ a iP) M[x 9l J + 2] . (6) 

Given Eq. © for 2 ff i(ic, Q 2 ), we choose the set 
{N m ax, ct, /?} to achieve optimal convergence of this series 
throughout the kinematic region constrained by the data. 
In practice, we find N max = 9, a = 3.0, and j3 = 0.5 to 
be sufficient. 



III. QCD ANALYSIS & PARAMETRIZATION 

A. Parameterization 

We consider a proton comprised of massless partons 
with helicity distributions q±(x, Q 2 ) which carry momen- 
tum fraction x with a characteristic scale Q. The differ- 
ence Sq(x,Q 2 ) = q + (x,Q 2 ) — q-(x,Q 2 ) measures how 
much the parton of flavor q "remembers" of the parent 
proton polarization. We will parameterize these polar- 
ized PDFs at initial scale Q 2 , = 4 GeV 2 using the follow- 
ing form: 

x 6q(x, Ql) = Afgiigx^ (1 - x) b * (I + c q x) , (7) 

where the polarized PDFs are determined by parameters 
{rj q , a q , b q , c q }, and the generic label q = {u v ,d v , q, g} de- 
notes the partonic flavors up-valence, down- valence, sea, 
and gluon, respectively. The normalization constants J\f q 

aW 1 + c %^) b( ^ + 1) ' (8) 

are chosen such that r\i are the first moments of 
5qt(x, Qq), rji = J Q dxSqt(x, Qq), where B(a, b) is the Eu- 
ler beta function. 

The total up and down PDFs are a sum of the valence 
plus sea distributions: Su = Su v + Sq and Sd = Sd v + Sq. 
We will assume an SU (3) flavor symmetry such that Sq = 
Su = Sd = Ss = Ss. While we could allow for an SU(3) 
symmetry violation term by introducing n such that Ss = 
5~s = nSq, as the strange PDF is poorly constrained the 
results would be insensitive to the specific choice of k. 

As seen from Eq. ([7]), each of four polarized par- 
ton densities q = {u Vl d v , q, g} contain four parameters 
{rj q , a q , b q , c q } which gives a total of 16 parameters that 
we must constrain. We now demonstrate that we can 
eliminate some of these parameters while maintaining 
sufficient flexibility to obtain a good fit. 
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1. First Moments of Su v and Sd v 

The parameters r] Uv and r]d v are the first moments of 
the 6u v and Sd v polarized valence quark densities; these 
quantities can be related to F and D as measured in 
neutron and hyperon /3-decays according to the relations 



a 3 = I dx Sq 3 = r) Uv - r)d v = F + D , (9) 
lo 



as 



f dx 5q 8 = r) Uv + f]d v =3F-D . (10) 
Jo 



where 03 and a% are non-singlet combinations of the first 
moments of the polarized parton densities corresponding 
to 

q 3 = (Su + Su) - {5d + Sd) , (11) 
q s = (5u + Su) + (5d + 8d)-2(5s + 5s) . (12) 

A reanalysis of F and D with updated /3-decay constants 
obtained (z3 F = 0.464 ± 0.008 and D = 0.806 ± 0.008. 
With these values we find: 



Vu v = +0.928 ± 0.014 , 
rj dv = -0.342 ±0.018 . 



(13) 
(14) 



We make use of r] Uv and 77^ to reduce the number of 
parameters by two. 



2. Gluon and Sea-Quarks 

We find the factor (l + c q x) in Eq. provides flexibil- 
ity to obtain a good description of the data, particularly 
for the valence distributions {«„, d v }. Thus we will make 
use of the c q coefficients for the the up- valence and down- 
valence distributions; in contrast, we are able to set the 
values for Cq and c g to zero (c§ = c g = 0) while main- 
taining a good fit and eliminating two free parameters. 
For the parameters {Cu v ,Cd v } we find the fit improves if 
we use non-zero values, but as these are relatively flat 
directions in %-space we shall fix the values as detailed 
in Table I. 

Separately, we find the b parameters control the large- a; 
behavior of the PDFs; thus, the sea-quark and gluon dis- 
tributions have large uncertainties in this region as they 
are dominated by the valence. To provide some guidance, 
we observe that for unpolarized parton densities in the 
large-x region, a ratio of b q /b g ~ 1.6 provides a good fit. 
Therefore we impose this ratio on the polarized bq and b g 
parameters to further reduce the free parameters. Addi- 
tionally, we are able to extract reasonable constraints on 
the a g and a g parameters; this is a benefit of the Jacobi 
polynomials. 

Having fixed {r] Uv , , Cq, c g } and the ratio b g /b g in 
preliminary minimization, we then set the parameters 



{bq,b g ,c u ^,Cdv\ as indicated in Table [TJ this gives us a 
total of 9 unknown parameters, in addition to (X s {Qq). 



B. DGLAP Evolution 

In the Jacobi polynomial approach the DGLAP evo- 
lution equations are solved in Mellin space. The Mellin 
transform of the parton densities q are defined analogous 
to that of Eq. ©: 

M[5q(x,Q 2 ),N] = Sq(N,Q 2 )= [ x N ^ Sq(x,Q 2 ) dx 

Jo 

N — l + a a 



' N + a q + b q / 
x B(N - I + a qi b q + 1) , (15) 

where q = {u v , d v ,q,g}, and B is the Euler beta function. 

In Mellin space, the twist-2 contributions to the polar- 
ized structure function gi(N,Q 2 ) can be represented in 
terms of the polarized parton densities and the coefficient 
functions ACf by: 

x [Sq(N,Q 2 ) + Sq(N,Q 2 )} 

+ ^2AC»6g(N,Q 2 )} . (16) 

Here, the sum runs over quark flavors {w, d, s}, and 
{<5g, 5q, Sg} are the polarized quark, anti-quark, and 
gluon distributions, respectively. 

The coefficient functions ACJj are the N-th moments 
of spin-dependent Wilson coefficients, and are given by 
El: 



ACf = |{- a(JV) + ( A(W+ (|-_^) 



AC, 



JV 



with Si {n) 



N - 1 
' N(N+ 1) 

En 1 
0=1 3 



1 



=1 F 



= (^) - V'(n + 1), ip(n) 



N(N 

h 7b, S 2 (n) = 
T'(n)/T(n) and 



ip'(n) = d 2 \nT(n)/dn 2 . 

In summary, we are able to express xg\ in terms of 9 
unknown parameters at an input scale of Qq = 4 GeV 2 . 
We now examine the fits to the spin structure functions 
to extract the polarized PDFs from the available data. 



IV. QCD FIT OF xgi(a;,Q 2 ) DATA 

Our analysis is performed using the QCD-PEGASUS 
program [7a|. We work at NLO in the QCD evolution 
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0.928 (fixed) 






-0.054 ± 0.029 




««„ 


0.535 ± 0.022 


8q 


dq 


0.474 ±0.121 




b Uv 


3.222 ± 0.085 




h 


9.310 (fixed) 




Cu v 


8.180 (fixed) 




Cg 









-0.342 (fixed) 




% 


0.224 ±0.118 


Sd v 


a d v 


0.530 ± 0.067 


5g 


"a 


2.833 ±0.528 




bd„ 


3.878 ± 0.451 




b g 


5.747 (fixed) 






4.789 (fixed) 




C 9 





a s (Qo) = 0.381 ±0.017 






X 2 /dof = 273.6/370 


= 0.74 



Table I: Final parameter values and their statistical errors in 
the MS-scheme at the input scale Qq — 4 GeV 2 . 



using Nf = 3 in the fixed-flavor number scheme with 
massless partonic flavors {u, d, s}. We take the renormal- 
ization and factorization scales to be equal (/ir = hf), 
and we compute the strong coupling a s (Q 2 ) at NLO us- 
ing a fourth order Runge-Kutta integration. Our initial 
parameterizations (Eq. [7|) are chosen to be invertible in 
N-space, and this makes our fitting procedure numeri- 
cally efficient. 

For the proton data we use EMC 0, HERMES (EEl, 
SMC [I, E143 @|, E155 [ll[ and COMPASSjll, for the 
neutron data we use £1420], HERMES (5l|ij~and E154 
d,0|, and for the deuteron data we use SMC (8|, E143 @, 
El 55 (T0| and HERMES [3. This data is summarized in 
Table HI 

We minimize the global X 2 [H, S [zi| : 



v 2 = V 

Aglobal / _j 



(17) 



where the sum n runs over the different experiments, w n 
is a weight factor for the n-th experiment, and x\ i s gi ven 
by: 



AA/"„ 



/ Kf exp _ theor 



exp 



(18) 
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Figure 1: The polarized parton distribution as function of x 
and for different values of Q 2 . 
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Figure 2: The polarized parton distribution at Q\ — 4 GeV 2 
as a function of x. Our fit is the solid curve. Also shown are 
the results of BB (dashed) Eel DSSV (dashed-dotted) [H], 
GRSV (long dashed-dotted) [28fl. and AAC (dashed-dashed- 
dotted) H|. 



Here, g^f ', Ag™ P j and g\ h l° r denote the experimental 
measurement, the experimental uncertainty (statistical 
and systematic combined in quadrature) and theoretical 
value for the i tli data point, respectively. AAf n is the ex- 
perimental normalization uncertainty and J\f n is an over- 
all normalization factor for the data of experiment n. We 
allow for a relative normalization shift J\f n between dif- 
ferent data sets within uncertainties AAf n quoted by the 
experiments. 

We minimize the above x 2 value with the 9 unknown 
parameters plus an undetermined a s (Qg). The values of 
these parameters are summarized in Table |TJ We find 
X 2 /d.o.f. = 273.6/370 which yields an acceptable fit to 
the experimental data. 



V. PDF AND STRUCTURE FUNCTION 
ANALYSIS 

We next present our polarized PDFs and perform com- 
parisons with other recent parameterizations @, [3lTl34| . 



A. Polarized PDFs 

Figure [1] displays our polarized PDFs for a selection of 
Q 2 values. The up- valence (xSu v ) and gluon (xSg) distri- 
butions are positive, while the down- valence (x8d v ) and 
sea (x6q) distributions are negative. We observe that the 
evolution shifts all the distributions to smaller values of 



Experiment 


rr-range 


Q^-rangelGeV 2 ] 


# of data points 


Mi 


E143 (p) 


0.031-0.749 


1.27-9.52 


28 


0.9998 


HERMES (p) 


0.028-0.66 


1.01-7.36 


39 


1.0006 


SMC (p) 


0.005-0.480 


1.30-58.0 


12 


0.9999 


EMC (p) 


0.015-0.466 


3.50-29.5 


10 


1.0094 


E155 


0.015-0.750 


1.22-34.72 


24 


1.0226 


HERMES06 (p) 


0.026-0.731 


1.12-14.29 


51 


0.9992 


COMPASSIO (p) 


0.005-0.568 


1.10-62.10 


15 


0.9920 


Proton 






179 




E143 (d) 


0.031-0.749 


1.27-9.52 


28 


0.9990 


E155 (d) 


0.015-0.750 


1.22-34.79 


24 


0.9998 


SMC (d) 


0.005-0.479 


1.30-54.80 


12 


0.9999 


HERMES06 (d) 


0.026-0.731 


1.12-14.29 


51 


0.9976 


Deuteron 






115 




E142 (n) 


0.035-0.466 


1.10-5.50 


8 


0.9991 


HERMES (n) 


0.033-0.464 


1.22-5.25 


9 


0.9999 


E154 (n) 


0.017-0.564 


1.20-15.00 


17 


0.9996 


HERMES06 (n) 


0.026-0.731 


1.12-14.29 


51 


1.0000 


Neutron 






85 




Total 






379 





Table II: Published data points with the measured x and Q 2 ranges, the number of data points (with a cut of Q 2 > 1.0 GeV 2 ), 
and the fitted normalization shifts Mi. 



x, and tends to flatten out the peak for increasing Q 2 . 
Figure [2] displays the extracted NLO polarized PDFs as 
compared with various parameterizations from the liter- 
ature mum!. 

Examining the xSu v and xdq distributions we see that 
most of the fits are in agreement, with the possible excep- 
tion of the DSSV [38| curves; for both distributions, the 
DSSV results approach zero more quickly than the other 
curves. For the x5d v distribution, all of the curves are 
comparable. The DSSV analysis employs results from 
semi-inclusive DIS (SI-DIS) data which can impose in- 
dividual constraints on individual quark flavor distribu- 
tions in the nucleon [38j. Finally, for the gluon distribu- 
tion, the DSSV results have a sign change in the region 
of x ~ 0.1 while the other fits are positive. Our result 
for gluon distribution is located between DSSV curve and 
the other fits (2Sl. l39l. l40l| . In particular, we find the gluon 
polarization vanished more quickly for small x values as 
compared with the other fits; we conjecture that using 
available asymmetry data in low x region may contribute 
to this difference. 



B. gi Structure Functions 

Figure [3] displays results for the polarized structure 
function xg\. For comparison, we display the results ob- 
tained by (Blumlein, Bottcher) BB [31|, (Gluck, Reya, 
Stratmann, Vogelsang) GRSV [HI, (Leader, Sidorov, 
Stamenov) LSS [H], (de Florian, Navarro, Sassot) DNS 
[HI and (Asymmetry Analysis Collaboration) AAC [3^ |. 
There is some spread in the analyses at low values of x; 
however, the data are generally well described within er- 
rors. As in the unpolarized case, the presence of scaling 
violations result a slope that varies with changing x val- 



ues; this is evident in Figure [3] where we observe the Q 2 
dependence of the structure function gi(x, Q 2 ). 

Given the polarized proton PDFs, we can use isospin 
symmetry to obtain the corresponding neutron structure 
functions. In Figure [4l we plot the neutron polarized 
structure function xg™. We also display the NLO QCD 
curves obtained by Ref. [44| in the polarized valon model 
(PVM). 

We can relate the deuteron structure function to that 
of the proton and neutron via: 

M[gt,N] = ± (l-ltj D ){M\£,N] + M\9?,X\) , 

(19) 

where ujjj = 0.05 ± 0.01 is the D-state wave probability 
for the deuteron [t3| - In Figure [S] we present our results 
for the structure functions xgi(x,Q 2 ), xgi(x,Q 2 ) and 
xgf(x, Q 2 ). and this compares favorably with the results 
of the BB [HI, GRSV [H, LSS [3, DNS [H and AAC 
[3^ | analyzes. 

The non-singlet spin structure function xg^ s (x, Q 2 ) is 
defined as [l2| 

xg? s (x,Q 2 ) ee xgl(x,Q 2 ) - xg"(x,Q 2 ) 

= 2[xgf(x,Q 2 )~ X f {x ; Q2) ]. (20) 

This is displayed in Figure [5J and we compare with the 
HERMES data [l2| for various Q 2 bins. In the second 
line of Eq. (|2"0|) we have related the structure function of 
the deuteron using isospin symmetry and the relation of 
Eq. (HD. 
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Figure 3: The polarized structure function g^ as function of 
Q 2 in intervals of x. The error bars shown are the statistical 
and systematic uncertainties added in quadrature. Our fit is 
the solid curve. The values of the shift a are given in paren- 
theses. Also shown are the results of BB (dashed) I31II. GRSV 
(dashed-dotted) HU, LSSfdashed-dotted-dotted) [3J], DNS 
(dashed-dashed-dotted) [33l | and AAC (long dashed-dotted) 
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Figure 4: The polarized structure function xg" as function of 
x and for a fixed value of Q 2 = 5 GeV 2 . The present fit is the 
solid curve. Also shown are the results of AK [44| (dashed) 
according to polarized valon model (PVM). 

~G? 2 GeV 2 3 GeV 2 5 GeV 2 10 GeV 2 



Att„ 
Ad v 
Aq 
Afl 

r? 
r; 
rf 



0.928864 
-0.342318 
-0.053400 
0.143610 
0.128291 
-0.050972 
0.035296 



0.928310 
-0.342114 
-0.053893 
0.191313 
0.131199 
-0.052416 
0.035965 



0.927794 
-0.341924 
-0.054379 
0.248845 
0.133822 
-0.053735 
0.036559 



0.927288 
-0.341738 
-0.054789 
0.323886 
0.136303 
-0.055000 
0.037115 



Table III: The first moments of polarized parton distributions, 
Au v , Ad v , Aq, Ag and polarized structure functions rf, r™, 
rf in NLO in the MS-scheme for some different values of Q 2 . 



D. First moment of g\ structure functions 



C. (?2 Structure Function 

We can now extract the structure function xg-i via the 
Wandzura-Wilczek relation [78l.l79l|: 



92{x, Q 2 ) = -gf(x, Q 2 ) + f %{(y, Q 2 ) . (21) 

Jx y 

This relation remains valid in the presence of target mass 
corrections. In Figure [7] we show our result for xg2 and 
we compare it with the experimental data from E143 [|J 
and SMC @|. 



We next use the polarized PDFs to compute the first 
moments, and compare with other recent analyzes. We 
can obtain the first moment of g± by 

T{(Q 2 )= [ dxg{{x,Q 2 )- (22) 
Jo 

The results of our fit are presented in Table IIIII for se- 
lected values of Q 2 , and these are compared with results 
from the literature in Table ITVl 

In the framework of QCD the spin of the proton can be 
expressed in terms of the first moment of the total quark 
and gluon helicity distributions and their orbital angular 
momentum, i.e. 

l - = 1a57 + AgP + LI , (23) 




Figure 5: The polarized structure function xgf, xg" and xgf as a function of x for selected values of Q 2 . The data are well 
described by the fit (solid curve). Also shown are the QCD NLO curves obtained by BB (dashed) [3lj . GRSV (dotted) [28| . 
LSS (dashed-dotted) [Hi, AAC (dashed-dotted-dotted) Q and DNS (dashed-dashed-dotted) [HJ. 





Model 


BB [40] 


GRSV [28] 


AAC [32] 


Ait i, 


0.928 


0.928 


0.9206 


0.9278 


Ad v 


-0.342 


-0.342 


-0.3409 


-0.3416 


Au 


0.874 


0.866 


0.8593 


0.8399 


Ad 


-0.396 


-0.404 


-0.4043 


-0.4295 


Aq 


-0.054 


-0.062 


-0.0625 


-0.0879 


A 3 


0.224 


0.462 


0.6828 


0.8076 



Table IV: Comparison of the first moments of the polarized 
parton densities in NLO in the MS-scheme at Q 2 — 4 GeV 2 
for different sets of recent parton parameterizations. The 
second column (Model) contains the first moments which is 
obtained from our new parametrization based on the Jacobi 
polynomials expansion method. The BB [13], GRSV and 
AAC results are also shown. 



typical value of Q 2 = 4 GeV 2 is around 0.355 in our 
analysis. We can also compare this value in NLO with 
other recent analysis. The reported value from the BB 
model [H is 0.569^116 AAC model [H] is 0.837 and 
the GRSV model ^ is 0.785, while the DSSV model 
[HI is approximately 0.1. Since the values of ^AS are 
comparable, we observe that the difference between the 
above reported values must come from different gluon 
distributions. 



E. Strong Coupling Constant 



where L p z is the total orbital angular momentum of all 
quarks and gluons. The contribution of ^AE + Ag for 



In this QCD analysis we extract a s (Q§) at NLO and 
obtain 
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HEREMES L12<Q <Z87 
MODEL Q 2 =1.5 GeV 2 




o.i 

X 



: • HERMES 6 32<Q<9.56 
O- 2 : MODEL Q 2 =8 GcV 2 



S> — 



0.1 
X 



HERMES 3.08<Q<5.6 
MODEL Q 2 =4J GeV 2 



S> 




0.1 
X 



: • HERMES 1L36<Q <1429 
O- 2 : MODEL Q 2 =13 GcV 2 

"a 



S> : 



0.1 
X 



Figure 6: The non-singlet polarized structure function xg^ s 
as function of x. 



ry (Mi) 


Order 


Reference 


Notes 

x i yj LCD 


0.1149 ± 0.0015 


NLO 




This analysis 


0-11321^™°° 


NLO 




40] 




„ „„.+o ooii) 
0-1134l ;™^ 


NNLO 




80] 




0.1141 ± 0.0036 


NLO 




44] 




0.1131 ± 0.0019 


NNLO 




63] 




0.1139 ± 0.0020 


l\ l\ T i\ t r s~~\ 

NNNLO 




66] 




1141 i u uuzu 


NNNLO 




80] 




0.1135 ± 0.0014 


NNLO 




811 


FFS 


0.1129 ± 0.0014 


NNLO 




811 


BSMN 


0.1124 ± 0.0020 


NNLO 




82] 


dynamic approach 


0.1158 ± 0.0035 


NNLO 




82] 


standard approach 


0.1171 ±0.0014 


NNLO 




83] 




0.1147 ±0.0012 


NNLO 




84] 




0.1145 ±0.0042 


NNLO 




85| 


(Preliminary) 


0.1184 ±0.0007 






86] 


World Average 



Table V: Comparison of a s (Mz) values from the literature. 
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Figure 7: The polarized structure function xg2 as function of 
x for Q 2 = 5 GeV 2 . 



n p 



a n 


= 5.650556817 


P 

% 


= 0.0148376 


K 


= 4.15388 


V 


= 0.986818274 


0? 


= -0.0189575 




= -4.75525 


c n 


= 0.064446823 


P 

Oa 


= 0.0121792 




= 2.68417 


d" 


= 0.807650292 


p 

a 3 


= -0.0040397 




= -0.800306 






P 

a 4 


= 0.000540845 




= 0.101095 



Table VI: Numerical coefficients for Eq sJ|28l29p of Af? He (y) 
and Afs H (y) obtained from Refs. [88H9(]| |. 



F. Nuclear Polarized Structure Functions 



Using the polarized PDF fit results, we examine the 
nucleon corrections factors for 3 He and 3 H. The polar- 
ized structure functions g l He and g x H can be composed 
from the polarized proton structure g^ and the polarized 
neutron structure gf as follows: 



a s {Ql) =0.381 ±0.017 . (24) 

Rescaling this to the Z boson mass scale we find 

a s {M%) = 0.1149 ± 0.0015 . (25) 

The error given in above does not include the relative 
systematics of the different classes of measurements. In 
Table [V] we provide a comparison of this value with other 
determinations from the literature computed at NLO 
and higher orders, including the current world average 
of a s {M%) = 0.1184 ± 0.0007. 




y 

Figure 8: The polarized light cone distribution function for 
the proton in the 3 He, based on the results of Ref. (88l - [90l |. 
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Figure 9: The polarized light cone distribution function for 



the neutron in the He, based on the results of Ref. [88 9 1 



0.5 



-0.5 



-1.5 



bV)l 



Q 2 =2.5 GeV 2 




NLO 






• E142 1.1<Q 2 <5.5 GeV 2 - 




A JLAB 2.71<Q 2 <4.83 




PVM 
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MODEL 



0.01 



0.1 



Figure 10: Analytical result for the polarized He structure 
function v.s. x for fixed Q 2 = 2.5 GeV 2 . The current fit is the 
solid curve. Also shown are the QCD NLO curves obtained by 
AK (dashed) according to polarized valon model (PVM) 
and BB (dashed-dotted) |3ll |. 




Figure 11: Analytical result for the polarized H structure 
function v.s. x for fixed Q 2 = 2.5 GeV 2 . The current fit is the 
solid curve. Also shown are the QCD NLO curves obtained by 
AK (dashed) according to polarized valon model (PVM) 
and BB (dashed-dotted) |3ll | for comparison. 



0.014[ 5 [(x,Q 2 )- 4.^(2-, Q 2 )] , (26) 



gl H (x,Q 2 ) = 2 f^A/^K^Q 2 ) 

J oo v v 

3 ^Afi H (y) 9 ^,Q*) 

+ O.OM^Q 2 )^^ 2 )] . (27) 

Here, Afi^^y) and Afi^ H (y) are the spin -dep endent nu- 
cleon light-cone momentum distributions |87l |88|. These 
functions parameterize the Fermi motion and the nucleon 
binding, and are readily calculated using the ground- 
state wave functions of 3 He and 3 H. Note that the last 
term in above equations is important only in the large- a; 
region. 

If we utilize isospin symmetry, we can equate Af% He (y) 
to Af? H (y), and also Af? He (y) to Af£ H (y); thus, we are 
left with only two independent functions Aff He (y) and 
Af3 He (y). Using the results of Refs. (sH-HO], we express 
these distributions as 



A/aV(y) 



0.5(l-d")(-i>"+j/)^ 

$*ja 



'(-b"+'») 2 



(28) 



(29) 



where U n (y) is a Chebyshev polynomials of the sec- 
ond type. The numerical coefficients of these equa- 
tions are presented in Table IVII We can then use 
Eqs. (|26l27p to obtain the polarized nucleon structure 
functions g 1 He (x,Q 2 ) and g l H {x,Q 2 ). 

To determine the g\ He and g± polarized structure 
functions we need the polarized light-cone distribution 
functions for proton and neutron in 3 He, i.e. Afi He and 
A/Jrj . In Figures [S] and IH1 we present our results using 
the parametrization of E qs. (I28I29P which is based on the 
numerical results of Ref. (901 • 

In Figures ITOl and fTTI we show our results for the g\ He 
and t?i H polarized structure function, and compare with 
BB [3l|, and the polarized valon model (PVM) 0]. For 
the gi He polarized structure function we see that our 
result coincides with the BB fit for x values down to 
~ 10~ 2 , and then falls off more quickly at very small x 
values. The polarized valon model (PVM), while still a 
reasonable fit to the data, lies below both of the other 
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fits. For the g x polarized structure function, our fit co- 
incides with the BB fit at both large and small x values, 
but dips below it (closer to the PVM) for intermediate x 
values. The differences between these curves come from 
the various data sets used, the constraints imposed, and 
the form of the parameterization. For example, in the 
AK fit [44] , only 257 experimental data points were used 
as the neutron data were not included; in contrast, the 
present analysis uses 379 points which does include the 
neutron data. Furthermore, the AK fit used 15 free pa- 
rameters while there are only 9 free parameters in the 
present analysis. These differences are reflected in the 
extractions of PPDFs, and a comparison of these differ- 
ent analyses may be indicative of the stability of the de- 
termined QCD parameters. 



mial method. Having extracted the polarized PDFs, we 
compute various nuclear structure functions ((71,(72) and 
Bjorken sum rule. In general, we find good agreement 
with the experimental data, and our results are in ac- 
cord with other determinations from the literature; col- 
lectively, this demonstrates progress of the field toward a 
detailed description of the spin structure of the nucleon. 

Having demonstrated the compatibility of the Jacobi 
polynomial method with other approaches in the litera- 
ture, this study can serve as a foundation for addressing 
issues of polarized scattering processes from a comple- 
mentary perspective. In particular, the Jacobi polyno- 
mial method offers the opportunity to examine efficien- 
cies of different methods, and this work is in progress. 



G. Bjorken Sum Rule 

We can also study the Bjorken sum rule [9l| which 
relates the difference of the first moments of the proton 
and neutron spin structure functions to the axial vector 
coupling constant of the neutron /3-decay, 

f [ g p 1 (x,Q 2 )-g?(x,Q 2 )}dx=lg A [l+0(^)] , (30) 

Jo 6 7T 

where g A = 1.2670 ± 0.0035 [zl, and the QCD radiative 
corrections are denoted as 0(— ). This sum rule can be 
generalized for the 3 He- 3 H system as follows: 

f\l H (x,Q 2 ) -gl He (x 7 Q 2 )]dx = \ ~g A [l+0(^)} , 
Jo b TT 

(31) 

where g A is the axial vector coupling constant of the Tri- 
ton /?-decay, with g A = 1.211 ± 0.002 Taking the 
ratio of the Eqs. $U\) and ([5T|). we find 

S^\9\ H M 2 )-9i He M 2 )W = 9a m 

fZ\£(x t Q*) - g?(x,Q 2 )]dx 9a- 

Given g A and g Al we compute the above ratio to be 0.956 
[88[. Note that the QCD radiative corrections are ex- 
pected to cancel exactly in above equation. Using the 
Bjorken sum rules of Eqs. (|30I31[) . we obtain the value 
0.924 for the ratio of Eq. 

VI. CONCLUSIONS 

We have presented a fit to the polarized lepton-DIS 
data on nuclei at NLO QCD using the Jacobi polyno- 
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Appendix: FORTRAN-code 

A FORTRAN package containing our gi(x,Q 2 ) polar- 
ized structure functions for {p, n, d, NS, 3 He, 3 H} and 
xg\(x,Q 2 \ as well as the polarized parton densities 
{u v d v ,g,q}. x5u v (x,Q 2 ), xSd v (x,Q 2 ), x8g(x 1 Q 2 ) and 
xSq(x,Q 2 ) at NLO in the MS-scheme can be found 
in http : //particles . ipm. ir/links/QCD. htm or ob- 
tained via e-mail from the authors. These functions are 
interpolated using cubic splines in Q 2 and a linear in- 
terpolation in log (Q 2 ).The package includes an example 
program to illustrate the use of the routines. 
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